Fast optimal frictionless atom cooling in harmonic traps 
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A method is proposed to cool down atoms in a harmonic trap without phase-space compression 
as in a perfectly slow adiabatic expansion, i.e., keeping the populations of the instantaneous initial 
and final levels invariant, but in a much shorter time. This may require that the harmonic trap 
becomes an expulsive parabolic potential in some time interval. The cooling times achieved are also 
shorter than previous minimal times using optimal-control bang-bang methods and real frequencies. 
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A fast adiabatic expansion in a short finite time looks 
like a contradiction in terms. An "adiabatic" process 
in quantum mechanics is a slow process where the sys- 
tem follows at all times the instantaneous eigenvalues and 
eigenstates of the time-dependent Hamiltonian. This is 
in a sense maximally efficient as the populations do not 
change, i.e. there is no heating or friction, but the price 
to pay is that the long times needed may render the pro- 
cess useless or even impossible to implement. Thus, a 
highly desirable goal is to prepare the same final states 
and energies of the adiabatic process in a given finite time 
t f , without necessarily following the instantaneous eigen- 
states along the way. We would also like the process to be 
robust with respect to arbitrary initial states. If fulfilled, 
this old goal [l| has important implications. In partic- 
ular, cooling without phase-space compression, which is 
all that is needed for many applications other than Bosc 
Einstein condensation, could be performed in fast cycles 
increasing, for example, the flux of cold atoms produced 
and the signal to noise ratio in an atomic clock 0], or 
in cold-atom pulsed beam experiments and related tech- 
nology This goal also includes as a particular case 
a long standing question in the fields of optimal control 
theory and finite time thermodynamics, namely, to opti- 
mize the passage between two thermal states of a system 
0> S S 0| ■ For time-dependent harmonic oscillators, 
minimal times have been established using "bang-bang" 
real-frequency processes believed up to now to be opti- 
mal @, in which the frequencies are changed suddenly 
at certain instants but kept constant otherwise. In this 
letter we shall describe a robust solution to the stated 
general goal for atoms trapped in a time-dependent har- 
monic oscillator which applies both to equilibrium and 
non-equilibrium states. In particular we describe cooling 
processes performed in a time interval smaller than the 
minimal time of the bang-bang methods considered so 
far. We shall for simplicity describe our method for states 
representing single atoms of mass m, but the same results 



are immediately applicable to A-body non-interacting 
fermions or to a Tonks-Girardeau gas [8| , and generaliza- 
tions will be relevant for other driving processes, such as 
cold atom launching or the transport of ultracold atoms 
with optical tweezers 

We consider an effectively one dimensional time de- 
pendent harmonic oscillator, H = p 2 /2m + muj(t) 2 q 2 /2, 
with an initial angular frequency u>(0) > at time t = 
and final frequency ujf = u>(tt) < w(0) at time tf. (This 
amounts to a temperature reduction by a factor ujf/oj(0) 
if the initial and final states are canonical.) The challenge 
is to find a trajectory u>(t) between these two values so 
that the populations of the oscillator levels n — 0,1,2... 
at tf are equal than the ones at t = 0. Our main tool to 
engineer u>(t) and the state dynamics will be the solution 
of the corresponding Schrodinger equation based on the 
existence of invariants of motion [J [13, 11, 13] of the form 
I{t) = l/2[{l/b 2 )q 2 muj 2 + ±tt 2 ], where tt = b(t)p - mbq 
plays the role of a momentum conjugate to q/b, the dots 
are derivatives with respect to time, and coq is in prin- 
ciple an arbitrary constant. The scaling, dimensionless 
function b — b(t) satisfies the subsidiary condition 



b + uj{t) 2 b = uji/b 3 , 



(1) 



an Ermakov equation where real solutions must be cho- 
sen to make I Hermitian [l3j |. luq is frequently rescaled 
to unity by a scale transformation of b [l|. Other con- 
venient choice is wo = w(0) as we shall see. I(t) has 
the structure of a harmonic oscillator Hamiltonian as 
well (as long as w 2 > 0), with time-dependent eigen- 
vectors \n(t)} and time-independent eigenvalues (n + 
l/2)Tiu!o. The general solution of the Schrodinger 
equation is a superposition of orthonormal "expanding 
modes" ip(t,x) = J2 n c n e M ™^ (x\n(t)) where a n {t) = 
— (n + l/2)wo /q dt' /b 2 , and the c„ are time independent 
amplitudes. For a single mode and uiq > 0, 
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The average position is zero and the standard devia- 
tion a — (J dxx 2 ^^ 2 ) 1 / 2 is proportional to b, a = 
b(n + 1/2) 1 / 2 / (muio/h) 1 / 2 , which underlines the physi- 
cal meaning of the scaling factor. 

A much studied case corresponds to the frequency 
scaling iM) = uo{Q)/b 2 with b = {At 2 + 2Bt + C) 1/2 
[ToL 11 . 03]. Substituting this in the subsidiary condi- 
tion gives ui 2 = co(0) 2 + AC — B 2 . For a trap with hard 
walls, the square-root-in-time scaling factor b (A = 0) has 
been shown to provide fast and efficient cooling [n| [n| . 
However, for harmonic traps, much more commonly re- 
alized in ultracold experiments, such time dependence 
leads to negative values of coq even for modest cooling ob- 
jectives. This makes Eq. ([2]) invalid and, moreover, lin- 
ear combinations of a continuum of non-square- integrable 
expanding modes would be needed to describe the evo- 
lution of any single eigenstate of the initial trap. This 
is of course only a drawback to calculate the dynamics, 
not to realize the expansion in the laboratory. Numerical 
results using other (adiabatic basis) methods [17( show 
that, even though the root-in-time scaling is singularly 
efficient for adiabatic following as discussed below, the 
cooling performance fails for very short expansion times 
tf. An alternative, successful strategy put forward here, 
inspired in inverse scatterin g te chniques for complex po- 
tential optimization 1|| 3, 13], is to leave u)(t) undeter- 
mined at first and impose properties on b and its deriva- 
tives at the boundaries, t = and tf, to assure: (a) 
that any eigenstate of H(0) evolves as a single expand- 
ing mode and that (b) this expanding mode becomes, up 
to a position-independent phase factor, equal to the cor- 
responding eigenstate of the Hamiltonian H(tf) of the 
final trap. This keeps the populations in the instanta- 
neous basis equal at the initial and final times. After 
b(t) and its derivatives are fixed at the boundaries, b(t) 
may be chosen as a real function satisfying the boundary 
conditions, for example as a polynomial or some other 
convenient functional form with enough free parameters. 
Once b{t) has been determined, the physical frequency 
uj(t) is obtained from the subsidiary equation (fTJ). 

Let us first discuss the conditions at t = 0. By choos- 
ing 6(0) = 1, 6(0) = 0, H(0) and 1(0) commute and 
have common eigenfunctions at that instant. We set 
ujq = uj(0) from now on so that 6(0) = must hold 
as well. These boundary conditions imply that any ini- 
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FIG. 1: (color online). Examples of ansatz for b: A simple 
polynomial ansatz (solid line, Eq. and an exponential of 

a polynomial (dashed line, exp 

E^ d 3 t J ')- ^(0) = 250 x 2tt 
Hz, u)(t f ) = 2.5 x 2tt Hz, 7 = 10. 



tial eigenstate of H(0), u n (0), will evolve according to 
the expanding mode @ for all later times. In general 
H(t) and I(t) will not commute for t > 0, so that the 
expanding mode ^ n (t) may have more than one com- 
ponent in the "adiabatic basis" of instantaneous eigen- 
states of H(t), {u n (t)},n = 0,1,2..., where u n (t) = 

' . e xp[-%-uj(t)x 2 ] H, 



muj{t) t 



At 



(2™n!) 1 /2 

time tf we want ^ n (tf) to be proportional, up to the 
global phase factor e l0tn ^ t f' , to the corresponding eigen- 
state of the final trap u n (tt). To this end we impose 
b(t f ) = 7 = [wo/w/] 1 / 2 , b{t f ) = 0, b{tf) = 0. From 
Eq. ([2|), one finds (H(tf)) n in terms of 6/ = b[tf) and 
6/ = db(t)/dt\t=t f ■ Since 6/ and 6/ can be set indepen- 
dently we can minimize the terms depending on them 
separately, and the global minimum is found to be pre- 
cisely at the adiabatic energy (n + l/2)TkVf, which cor- 
responds to our boundary conditions. Any other choice 
would necessarily produce "frictional heating" . 
Substituting the simple polynomial ansatz 



b(t) 



3=0 



(3) 



into the six boundary conditions gives six equations 
that can be solved to provide the coefficients, 6(t) = 
6 (7 - 1) s 5 - 15 (7 - 1) s 4 + 10 (7 - 1) s 3 + 1, where s = 
t/tf, see Fig. [T] At initial and final times and tf, 
oj(t) = ujo/b(t) 2 , but this relation does not hold in gen- 
eral for an arbitrary intermediate time. 

The above mentioned six conditions leave time depen- 
dent phases e i<Xn ^ t > of no relevance regarding the popu- 
lation of the nth level. In particular stationary density 
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operators with respect to H(0) (e.g. a canonical state, 
or a pure state \u n (0)) (u n (0)|) are mapped onto the cor- 
responding stationary states of H{tf) with the phases 
canceled. In other cases the phases remain but the pop- 
ulations are preserved. Note that e lan ^ l \ see Eq. @, 
is the phase factor that the initial state u n (Q) would ac- 
quire in a virtual adiabatic process in which the adiabatic 
(instantaneous) energy had the form (n + l/2)fwj /b 2 . 
Phase control may as well be imposed by adding integral 
conditions such as r(tf) = J Q f dty-yz = ^t', where t' 
is some desired time. This of course requires an ansatz 
more complicated than Eq. ([3]), such as a polynomial of 
higher degree. 

Numerical examples of frequencies tu(t) and energies 
(H(t)) of fast adiabatic-like expansions are provided in 
Figs. 2-4 using the b shown in Fig. [T]for ujo ~ 250 x 2tt Hz 
and ujf = 2.5 x 2n Hz ( 7 = 10) . These values can be found 
in actual experiments [21| . We could formally study sub- 
hertz frequencies lo j but they would render the trap very 
sensitive to low- frequency acoustic noise [22]. Compare 
first the finite times considered (from 2 to 25 ms) with the 
times necessary for actual adiabatic following during the 
whole interval < t < tf. The adiabaticity condition for 
the harmonic oscillator becomes \^/2cu / (8uj 2 )\ <C 1. For 
a linear ramp, co(t) — ► luq + (ujf — LJo)t/tf, this implies a 
very long time, tf 3> 1.1 s. In fact it would be necessary 
to expand the trap for 6 s to achieve a 1% relative error 
in the final energy of the ground state. A much more ef- 
ficient strategy is to distribute uj/lo uniformly along the 



trajectory, i.e., lo/uj 



c, c 



being constant. Thus the 



trap expansion speed decreases with the splitting. Solv- 
ing this differential equation and imposing ujf = Lo(tf) we 
get uj{t) = wo/[l _ (j^f — too)t/(tfujf)}. This corresponds 
to the case A = 0, 2B — — {uif — ojq)/ '(tfUif),C — 1 
(i.e., a square-root-in-time scaling factor), and implies 
tf ^> 11ms from the adiabaticity condition. A 1% error 
level for the ground state energy is achieved after 45 ms. 

A prominent feature of the trajectories, see Fig. 2b, is 
that oj(t) 2 may be negative during some time interval in 
which the potential becomes an expulsive parabola. This 
is physically feasible and has been realized experimen- 
tally with an offset magnetic field that overcomes the op- 
tical dipole well in the axial direction of elongated cigar- 
shaped optical traps [23| . In general the (imaginary) fre- 
quency of the repulsive region increases for shorter cool- 
ing times as shown in Fig. [2j3. 

The appearance of transient energies below the final 
one, see e.g. the solid line in Fig. [3^ near t/tf = 0.15, 
may be misleading. It is a consequence of the repulsive 
regime and should not be interpreted as useful cooling 
in a time shorter than tf. Since the "trap" is actually 
a repeller the kinetic energy would grow without bound 
if the potential were kept frozen at the time when the 
energy is minimal. Similarly, if the potential were sud- 
denly changed into its final form, V(tf), the total energy 
would be higher than the adiabatic energy, i.e., the one 
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FIG. 2: (color online), (a) The average energies of the ground 
state expanding mode for different final times tf. tf =25 ms 
(solid), tf — 15ms (dashed), tf = 10ms (dotted), and 
tf — 6 ms (dash-dotted). Other parameters as in Fig. [T] (poly- 
nomial 6) (b) The corresponding squared frequency u>(t) 2 . 
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FIG. 3: (color online). Cooling in tf = 2ms. (a) Average 
energies of the ground expanding mode for b taken as a poly- 
nomial (solid), as an exponential of a polynomial (dashed), 
and for a piecewise constant frequency "bang-bang" process 
(dot-dashed) with uj\ = i0.9tJo and u>2 = loo- Other parame- 
ters as in Fig. [1] (b) The squared frequencies uj(t) 2 . 



for a population-preserving process. 

Fig. [SI illustrates that a given cooling objective may 
be attained in less time than the minimal time required 
by real-frequency bang-bang trajectories, optimal among 
real- frequency trajectories [6J. For the three-jump "tra- 
jectory" 0, 



uj (t = 0) 

Ul (0 < t < Tl) 

lo 2 (n < t < ti + t 2 ) 

Uf (t = tf = n + t 2 ) 



(4) 



the smaller u>i and the larger L02 are, the faster the cooling 
is. Thus the fastest process to reach the target state 
corresponds to the limit of toi — > and u>2 — > 00 [6| with 



y/1 - Uf/UJQ 



(5) 



These results are based on optimal control theory, initial 
and final thermal states, and the constraint wi,2 > 0. 
Clearly, relaxing the positivity condition for the interme- 
diate frequencies, makes faster processes with tf < tj ,m 
possible, which, moreover, involve only finite frequencies. 
Since Eq. has been used to justify a finite time version 



of the third principle (if LOf — > 0, 
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FIG. 4: (color online) . Average energies for expanding modes 
n = 1 (solid), n = 2 (dashed), and n = 3 (dotted), (a) 
t/ = 2 ms; (b) tf = 25 ms. Other parameters as in Fig. [T] 
(polynomial b). 



a revision of these conclusions. A bang-bang example is 
shown in Fig. 3 (dot-dashed lines), for tf = 2ms, much 
shorter than the time tff 1 ™ » 6 ms corresponding to the 
initial and final frequencies chosen. io\ = iuij is imagi- 
nary, and the corresponding b\{t) solving the Ermakov 
equation with initial conditions &i(0) = 1 and &i(0) = 1, 



takes the form bi(t) = [1 



sinh 2 (w/i)] 1/2 . In the 



second segment we assume u> 2 real which gives, with the 
final conditions b 2 (tf) = 7 and b 2 (tf) = 0, b 2 (t) = 



r 



sin 2 [0J2 (t-tf)]} 1 / 2 . The matching 

conditions b\ (ri) = &2(t"i) and b\ (t}) = 62 (ri) are then 
solved for n and i/, see Fig. 3 and its caption for details. 
Of course the discontinuous jumps in this type of trajec- 
tory call into question its realizability. Fig. 3 also shows 
two smooth trajectories for tf = 2 ms corresponding to 
two different ansatz for b (polynomial and exponential 
of a polynomial). The resulting rates of change of ui(t) 
are feasible with present technology. Indeed the intensity 
of a dipole trap (the frequency scales as the square root 
of the intensity) can be changed by three or four orders 
of magnitude in 100 ns using acousto-optics or electro- 
optics modulators. To monitor the sign of the square fre- 
quencies, one can superimpose two dipole beams locked 
respectively on the blue and red side of the line. By al- 
ternating them rapidly with a control of their relative 
intensity, one can shape the square frequencies and their 
signs at will. Alternatively, one can combine magnetic 
and dipole traps. 

Finally, while the initial (and final) state in Figs. 2 
and 3 is the ground state, Fig. 4 illustrates that the same 
uj(t) trajectories work as well for arbitrary excited states. 
This also means that fast frictionless cooling is directly 
applicable to arbitrary superpositions or mixed states, as 
well as to simple many body systems such as a polarized 
Fermi gas, or its (symmetrized) bosonic counterpart, the 
Tonks- Girardeau gas. Similar techniques may also be 
applicable for weakly interacting bosonic systems, the 
control of soliton dynamics of Bose-Einstein condensates 
and to manipulate the transport of ultracold 
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tionless cooling in a harmonic trap based on a modulation 
of the trap frequency u>(t) that includes in general time 
intervals in which the potential becomes expulsive. It is 
applicable to arbitrary initial states, not necessarily in 
equilibrium. The algorithm to design uj(t) depends on 
flexible ansatz functions that can be modified or made 
more and more complex, by adding parameters, to sat- 
isfy further requirements. This could be used, e.g., to 
minimize the maximal frequencies along the trajectory. 
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atoms [9J. 

We have in summary described a method for fast, fric- 
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